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HIGHLIGHTS 


•  Li-ion  cell  under  discharge  generates  thermal  energy  through  five  different  modes. 

•  Three  of  the  five  modes  are  the  electrolyte,  anode,  and  cathode  resistances. 

•  Two  of  the  five  modes  are  the  entropy  change  in  the  cathode  and  the  anode. 

•  The  five  are  dependent  on  the  rate  of  discharge  and  the  environmental  temperature. 

•  Measurement  and  quantification  of  heat  generated  by  each  mode  is  demonstrated. 
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A  lithium-ion  cell  under  discharge  generates  thermal  energy  (Q)  through  five  different  internal  pa¬ 
rameters  or  modes:  the  electrolyte  resistance  (Rs),  anode  resistance  (Ra),  cathode  resistance  (Rc).  and 
entropy  changes  in  the  cathode  (ASC),  and  the  anode  (ASa).  This  work  demonstrates  a  set  of  tools  to 
measure/quantify  the  heat  generated  by  each  parameter  separately  during  discharge.  These  five  sources 
are  not  dependent  upon  each  other;  they  are  dependent  on  the  state  of  charge  and  the  environmental 
temperature  (Tenv).  The  Q  generated  by  each  mode  varies  with  degree  of  discharge  and  Tenv.  Rs  generates 
most  of  the  Q  in  the  -10  °C  to  40  °C  range;  Rc  becomes  significant  at  Tenv  <20  °C.  Constant  current 
discharge  does  not  cause  a  monotonic  increase  in  anode  and  cathode  temperatures  (Ta  and  Tc),  due  to  the 
direction  of  change  in  A Sc  and  A Sa.  Negative  change  in  ASa  for  the  carbon  anode  cools  it,  causing  the  Ta  to 
level  off  and  even  decrease  with  increased  discharge.  A Sc  for  lithium  manganese  oxide  cathode  is  positive 
at  some  SoC  and  negative  at  others,  preventing  a  monotonic  increase  in  Tc.  Measuring  the  five  Qs 
separately  opens  the  opportunity  to  study  thermal-runaway  from  the  perspective  of  the  anode,  cathode 
and  electrolyte. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Lithium-ion  cells  are  among  the  optimal  electrical  energy  stor¬ 
age  and  power  delivery  devices.  During  each  charge— discharge 
cycle,  the  cell  converts  electrical  energy  into  chemical  and  thermal 
energies;  the  cell  is  only  capable  of  storing  the  chemical  energy.  The 
discharge  reverses  the  process,  converting  the  chemical  energy  into 
electrical  and  thermal  energies.  The  thermal  energy  generated  by 
the  cell  during  charging  and  discharging  is  permanently  lost,  never 
stored  or  recovered  as  electrical  or  chemical  energies,  and  accounts 
for  the  energy  loss  experienced  during  every  charge-discharge 
cycle.  Thermal  energy  generated  inside  the  cell  increases  the 
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internal  temperature  of  the  cell  and  under  extreme  conditions  can 
initiate  other  exothermic  chemical  reactions  causing  the  cell  to 
vent,  and  sometimes  transitioning  the  cell  to  a  thermal  runaway 
condition.  Venting  and  thermal  runaway  are  major  safety  issues 
related  to  Li-ion  cells,  and  their  causes  and  ramifications  can  be 
found  in  recent  reviews  and  news  articles  [1,23], 

The  most  common  and  accepted  modes  of  thermal  energy  (Q) 
generation  during  charging  and  discharging  of  Li-ion  cells  are  en¬ 
tropy-  and  internal  resistance-related,  and  researchers  use  indirect 
measurements  and  thermal  models  to  estimate  the  magnitude  of  Q 
[4-6],  The  more  direct  measurements  of  Q  are  made  using  calo¬ 
rimeters  [7—9],  However,  such  measurements  only  yield  the  sum  of 
all  the  heats,  not  the  individual  components  of  the  heat  generated 
inside  the  cell. 

It  has  been  commonly  recognized  that  the  heat  is  generated 
through  five  different  modes  by  five  independent  components 
inside  the  cell  [5—7],  They  are:  i)  the  electrolyte  resistance  (Rs); 
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ii)  the  anode  resistance  (Ra);  iii)  the  cathode  resistance  (Rc);  iv) 
the  entropy  change  in  the  cathode  (ASC);  and  v)  the  entropy 
change  in  the  anode  (ASa).  The  rate  of  heat  generated  by  each  one 
of  these  five  modes  is  dependent  upon  the  temperature  of  the 
cell,  and  the  amplitude  and  direction  of  the  current.  The  ability  to 
measure  the  individual  contribution  by  each  mode  to  Q,  espe¬ 
cially  during  the  course  of  charging  and  discharging,  yields 
invaluable  information  about  the  thermal  behavior  of  each  mode. 
It  is,  perhaps,  the  most  successful  way  to  track  all  the  vital  points 
inside  the  cell,  and  provide  a  trajectory  for  any  impending  errant 
behavior  by  the  cell. 

The  major  focus  of  the  work  described  below  is  to  quantify  the 
contributions  by  the  individual  modes  to  the  thermal  energy 
outputted  during  discharge,  and  to  reconcile  those  individual 
contributions  with  the  overall  thermal  behavior  of  the  cell  during 
discharge.  In  addition,  the  work  also  demonstrates:  i)  the  instan- 
taneity  of  the  impact  of  the  changes  in  the  discharge  current  on  Ta 
and  Tc;  ii)  ultraslow  responses  of  the  surface  temperature  (Tsurf)  to 
the  changes  in  the  current;  and  iii)  inhomogeneity  in  the  temper¬ 
ature  within  the  cell,  namely,  differences  between  the  Ta  and  Tc 


during  discharge;  and  iv)  presence  of  correlation  between  Ta  and  Tc 
with  Tenv,  and  absence  of  correlation  between  the  TSUrf  and  renv 

1.1.  Highlights  of  the  techniques  used  in  online  monitoring  of  heat 
generation 

Two  different  techniques  are  available  to  estimate  the  heat  en¬ 
ergy  generated  by  each  of  the  five  modes.  For  estimating  the 
entropy-related  thermal  energy,  we  used  a  recently-acquired 
technique  that  allows  one  to  measure  Ta  and  Tc  in  situ,  while  the 
cell  is  under  charge  and  discharge  [10],  By  measuring,  for  example, 
Tc,  one  could  compute  Qasc,  the  heat  associated  with  the  entropy 
change  in  the  cathode  (A Sc);  similarly,  Qas2  by  measuring  Ta.  On  the 
other  hand,  for  estimating  the  Joule  heating,  one  would  need  the  R3, 
Rc,  and  Rs  values  measured  during  discharge. 

The  Ra,  Rc,  and  Rs  values  can  be  obtained  using  conventional 
electrochemical  impedance  spectroscopy  (EIS).  The  technique  to 
measure  Ta  and  Tc  is  also  based  on  impedance.  A  measure  of  the 
phase  shift  (0)  between  an  applied  voltage  (or  current)  and  the 
resultant  current  (or  voltage)  at  a  single  frequency  is  a  function  of 
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the  temperature;  for  most  Li-ion  cells  tested  (capacity:  l-Ah-55- 
Ah),  the  ideal  frequency  for  measuring  Tc  is  in  the  10-  to  30-Hz 
range,  and  for  Ta  in  the  40-  to  100-Hz  range  [11 1, 

Note  that  the  capacitive  and  resistive  properties  of  the  anode 
and  the  cathode  are  quite  different  from  each  other,  such  that  0  is 
responsive  to  Ta  in  the  high-frequency  range  and  to  Tc  in  the  low- 
frequency  range.  Most  importantly,  the  frequency  range  where  0 
is  only  sensitive  to  Ta  is  well  separated  from  the  range  where  0  is 
only  sensitive  to  Tc ;  and  over  these  two  ranges,  the  phase  values  are 
insensitive  to  the  cell’s  state-of-charge  (SoC).  More  specifically,  at 
any  one  frequency  in  the  40-  to  100-Hz  range,  0  is  responsive  to  Ta; 
and  in  the  10-  to  30-Hz  range,  it  is  responsive  to  Tc  [11—13], 

The  impedance-based  approach  to  Ta  and  Tc  measurements 
provides  several  unique  advantages  over  the  fiber  optic-based 
internal  temperature  sensors  that  are  currently  under  develop¬ 
ment  [14],  The  fiber  optic  sensor  is  invasive.  If  it  is  too  thick 
(several  microns),  it  will  seriously  interfere  with  the  current 
density  distribution  between  the  anode  and  the  cathode,  and  if  it 
is  too  thin,  it  will  be  too  quickly  etched  by  the  multiple  chemicals 


(oxidizers  and  reducers)  that  are  present  inside  the  Li-ion  cell.  By 
contrast,  the  impedance-based  phase  measurement  is  truly 
noninvasive:  it  is  an  electrical  technique,  and  it  uses  electrical 
wires  to  connect  to  the  positive  and  negative  terminals  of  the  cell 
to  monitor  Ta  and  Tc.  Furthermore,  the  technique  facilitates 
measuring  Ta  and  Tc  separately  therefore  the  thermal  properties  of 
the  anode  and  cathode  can  be  assessed  separately  or  independent 
of  each  other. 

2.  Methods 

Measurements  of  heat  generation  during  constant-current 
discharge  at  C-rate  (4.4-A  rate)  were  conducted  on  freshly  pro¬ 
cured  4.4  Ah  Boston  Power  Swing  4400  cells.  Electrical  contacts 
to  the  cell  were  made  through  nickel  tabs  welded  to  the  cell’s 
positive  and  negative  terminals.  The  anode  in  the  cell  is 
composed  of  graphite;  the  cathode  is  lithium  manganese  oxide. 
The  freshly  obtained  cells  were  first  subject  to  two  or  three 
charge-discharge  cycles.  Next,  the  cells’  capacities  at  different 
temperatures  were  determined  by  cycling  them  at  the 
manufacturer-specified  charge— discharge  rates  using  a  pro¬ 
grammable  battery  tester  (Arbin  Instruments  Model  BT2000). 
During  the  tests,  the  cells  were  housed  inside  an  environmental 
chamber  (Tenney  Model:  TJR);  air  circulation  inside  the  chamber 
facilitated  active  heat  transfer  between  the  cell  surface  and  the 
surrounding.  A  surface-mounted  K-type  thermocouple  moni¬ 
tored  the  cell’s  surface  temperature,  Tsurf.  Prior  to  each  discharge, 
the  cells  were  first  discharged  to  2.7  V  cut-off  and  charged  to  full 
capacity  at  room  temperature  (23.5  °C),  and  next  kept  up  to  5  h 
in  an  environmental  thermal  chamber  at  the  designated  Tenv  in 
order  to  ensure  equilibration  between  the  cell’s  internal  tem¬ 
perature  Tenv.  Finally,  after  the  discharge  to  the  2.7  V  cut-off  value 
at  each  Tenv,  the  cell  temperature  was  equilibrated  with  the  room 
temperature  over  12  h,  and  re-discharged  at  C/4  rate  from  its 
open  circuit  voltage  to  a  2.7  V  cut-off;  this  procedure  yielded  the 
residual  electrical  energy  (Ed2)  that  was  left  un-discharged  dur¬ 
ing  the  first  discharge  at  the  Tenv.  The  heat  energy  generated 
during  the  second  discharge  (for  estimating  Ed2)  was  assumed 
negligible. 

During  the  course  of  charge  and  discharge,  the  cell  impedance 
(EIS)  and  phase  shifts  were  measured  using  a  Solartron  Model 
Electrochemical  Interface  (potentiostat/galvanostat)  Model  SI  1287, 
and  Frequency  Response  Analyzer  Model  SI  1250.  The  frequency  (f) 
range  for  the  EIS  measurement  was  set  to  span  from  0.8  Hz  to 
10  kHz,  for  a  total  of  45  individual  frequencies.  The  amplitude  of  the 
perturbation  current,  iac  was  set  at  50  mA  rms;  this  perturbation 
current  caused  no  more  than  a  1-mV  polarization  to  the  cell 
voltage,  Vac.  Ta  and  Tc,  as  well  as  Rs,  Ra  and  Rc  values  were  derived 
from  the  EIS  data. 

3.  Results  and  discussions 

The  discussions  below  are  centered  on  the  measurement  of 
thermal  energy  generated  during  discharge  by  five  independent 
modes:  Rs,  Ra,  Rc,  A Sc  and  ASa;  and  their  dependence  on  Tenv  The 
ability  to  measure  Ta  and  Tc  during  discharge  allowed  us  to  estimate 
Qasc  and  Qa s„-  Similarly,  the  ability  to  monitor  Rs,  Ra  and  Rc  during 
the  course  of  discharge  enabled  us  to  estimate  the  individual 
contributions  by  the  three  modes  to  Joule  heating.  The  impedance 
technique,  in  conjunction  with  the  use  of  electrical  equivalent  cir¬ 
cuits  to  measure  Rs,  Ra  and  Rc,  is  a  classical  method  to  probe  Li-ion 
cells,  and  is  amply  described  in  textbooks  [15],  The  technique  to 
measure  Ta  and  Tc  is  relatively  new  therefore  important  details  of 
the  technique  along  with  the  results  are  described  in  the  next 
subsection. 
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Fig.  3.  Internal  temperature  data  collected  during  the  discharge  of  the  4.4-Ah  cell  at  4.4-A  rate  in  six  different  Tenv  from  -10  °C  to  40  °C  as  indicated  in  (a)-(f).  The  dots  are  the 
experimental  Tc  and  Ta  values  measured  in  intervals  of  129.4  s.  The  dotted  lines  are  polynomial  fits  to  the  experimental  data. 


3.1.  Phase  shift  versus  Ta  and  Tc 

The  technique  used  in  measuring  Ta  and  Tc  is  based  on  electrical 
phase  (0)  measurement.  The  unique  relationships  between  0  and 
Ta;  and  <f>  and  Tc  are  demonstrated  in  Fig.  1  through  the  phase  shift 
versus  frequency  (0  vs.  f)  data  collected  at  20%,  40%,  60%  and  90% 
SoC  at  different  renv  in  the  -10  °C  to  40  °C  range.  Note  that  while 
collecting  the  data  shown  in  the  figure,  the  cell  was  not  under 
charge  or  discharge  in  order  to  avoid  current-induced  heat  gener¬ 
ation  inside  the  cell.  We  also  assumed  that  the  sub-microwatt  heat 
generated  by  the  50-mA  ac  current  used  in  perturbing  the  cell 
during  the  < p  measurement  is  too  small  to  affect  Ta  and  Tc.  Prior  to 
each  0  measurement,  the  cell  was  allowed  to  equilibrate  with  Tenv 
for  two  to  five  hours. 

At  each  temperature,  in  the  10-Hz— 700-Hz  range,  the  values  of 
0  are  found  to  be  virtually  independent  of  SoC  and  strongly 
dependent  on  Tenv ;  at/<  10-Hz,  the  0  values  are  dependent  on  SoC 
therefore  the  phase  data  at/<  10-Hz  should  not  be  used  to  measure 
Ta  or  Tc.  Note  that  the  data  shown  in  the  figure  were  collected  by 
first  setting  the  cell  at  a  given  SoC,  then  allowing  its  internal 


temperature  to  equilibrate  with  Tenv  in  a  thermal  chamber  for 
several  hours,  and  then  measuring  the  0  values;  most  importantly, 
after  setting  the  SoC  at  a  fixed  value,  for  the  entire  duration  of 
temperature  equilibration  and  0  measurements,  the  cell  was 
guarded  against  charging  or  discharging.  To  minimize  self¬ 
discharge  loss  of  stored  energy,  all  tests  were  completed  within  2 
days  of  charging,  Next,  an  arithmetic  mean  of  0  at  multiple  SoC  in 
the  5%— 95%  range  provided  a  single  data-point  in  a  calibration 
curve  of  “0  versus  Ta”  (for  example,  at/—  100  Hz)  and  “0  versus  Tc” 
(for  example,  at  /=  20  Hz).  Repeating  such  measurements  over  a 
wide  range  of  temperatures  range  constitutes  a  complete  0  versus 
Ta  and  0  versus  Tc  calibration  curves.  More  details  of  this  procedure 
to  obtain  0  versus  Ta  calibration  curve  are  described  in  Ref.  [11  ].  The 
0  vs.  Ta  and  0  vs.  Tc  calibration  curves  made  on  a  new  Swing  4400 
cell  are  shown  in  Fig.  2.  Note  that  one  set  each  of  calibration  for  Ta 
and  Tc  is  sufficient  for  given  cell  model  to  measure  its  internal 
temperatures,  Ta  and  Tc  under  static  (no  charge  or  discharge)  and 
dynamic  (during  charging  and  discharging)  conditions. 

Fig.  3  shows  the  Ta  and  Tc  values  estimated  using  0  measure¬ 
ments,  in  intervals  of  129.4  s  while  discharging  the  cell  at  4.4-A 
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Fig.  4.  Impedance  of  the  4.4-Ah  Li-ion  cell  at  various  Tenv  as  indicated.  For  clarity,  only 
one  set  of  impedance  at  one  SoC  (60%)  at  each  of  the  six  Tenv  are  shown  in  two  figures 
(a)  and  (b).  Each  point  in  the  graph  represents  impedance  data  at  one  of  about  45 
separate  frequencies, /between  0.82  and  1000  Hz  (or  150  Hz);  the  end  frequencies  are 
shown  in  the  figure.  The  dotted  lines  are  mathematical  fits  to  the  experimental  data, 
obtained  using  the  electrical  equivalent  circuit  shown  in  the  inset  in  (a).  The 
approximate  data  points  in  each  graph,  representing  the  cathode  and  the  anode  im¬ 
pedances,  are  shown  using  double  brackets. 

rate;  note  that  prior  to  each  discharge,  in  order  to  ensure  accuracy 
and  reproducibility  of  the  impedance  and  0  measurements,  the  test 
cell  was  first  fully  charged,  and  then  its  internal  temperature  was 
equilibrated  with  Tenv  for  several  hours  in  an  environmental  ther¬ 
mal  chamber. 

3.2.  Ra,  Rc  and  Rs  estimations  from  impedance  measurement 

During  discharging,  a  cell’s  SoC  and  the  temperatures  of  its  in¬ 
ternal  components  (anode,  cathode  and  the  electrolyte)  do  change. 
Our  goal  is  to  estimate  the  individual  thermal  energy  generated 
through  Joule  heating  by  Ra,  Rc  and  Rs  while  the  cell  is  under 
discharge.  A  feasible  path  to  the  estimation  is  through  a  measure  of 
the  cell’s  internal  impedance  at  multiple  frequencies  (in  the  1-Hz— 
10-kHz  range),  under  the  dynamic  condition  of  discharge.  Note  that 
at  the  discharge  rate  of  4.4-A,  a  fully-charged  4.4-Ah  cell  will 
discharge  in  roughly  1  h.  During  the  course  of  discharge  at  each  Tenv, 
we  made  about  30  sets  of  impedance  measurements  at  the  rate  of 
about  129.4  s  per  set.  One  such  set  at  each  Tenv,  for  a  total  six  sets, 
are  shown  in  Fig.  4.  Also  shown  in  the  figure  are  the  mathematical 
fits  to  the  data  at  each  Tenv,  and  the  electrical  equivalent  circuit  used 


in  arriving  at  those  fits.  The  fit  yields  Ra,  Rc  and  Rs,  and  by  repeating 
the  fit  over  all  the  sets  at  Tenv,  (roughly  a  total  of  180  sets  of 
impedance  data),  we  obtained  about  500  values  of  Ra,  Rc  and  Rs  at 
the  six  different  Tenv;  the  results  are  shown  in  Fig.  5. 

Note  that  Ra  has  a  weak  dependence  on  SoC  [16];  Rc  has  a 
stronger  dependence  on  SoC,  while  Rs  is  also  weakly  dependent  on 
SoC,  as  may  be  seen  from  Figs.  1  and  2  in  Ref.  [11].  Due  to  the 
varying  dependence  of  Ra,  Rc  and  Rs  on  SoC,  in  Fig.  5  and  in  several 
subsequent  figures  that  involve  the  three  resistances,  we  chose  to 
represent  their  dependence  on  temperature  as  a  function  of 
discharge  time. 

3.3.  ASa  and  ASC  estimations  from  literature  data 

Entropy- related  thermal  energy  (Qas)  is  obtainable  through  the 
relationship  Qas  =  —TAS.  Thus,  to  estimate  Qas  due  to  the  anode  and 
cathode,  we  need  two  parameters  for  each  electrode:  T  and  AS. 
Neither  T  nor  AS  remain  constant  during  the  course  of  discharge 
therefore  should  be  obtained  throughout  the  course  of  discharge. 
The  procedure  to  obtain  Ta  and  Tc  values  online,  i.e.  while  dis¬ 
charging  the  cell,  was  described  in  Section  3.1.  The  technique 
available  to  measure  AS  is  much  slower  than  impedance  mea¬ 
surement  therefore  monitoring  it  during  discharge  is  virtually 
infeasible;  details  about  estimating  AS  at  different  SoC  (under  static 
conditions  only:  no  charging  or  discharging)  are  described  in  the 
published  literature  [5,17—19],  Fortunately,  AS  is  not  dependent  on 
discharge  rate  or  renv;  it  is  only  dependent  on  SoC.  The  published 
literature  is  abundant  with  ASa  for  various  types  of  carbon,  the 
anode  material,  and  ASC  for  different  types  of  cathode  material. 

Part  of  the  cathode  material  in  the  Swing  4400  cell  is  lithium 
manganese  oxide  (LMO);  a  closely  matching  material  is  lithium  in 
spinel  manganese  oxide  for  which  the  ASC  vs.  SoC  due  to  the 
intercalation  of  lithium  was  estimated  by  Yazami  et  al.  [18]  The 
anode  material  in  the  Swing  4400  cell  is  graphite;  we  assumed  the 
ASa  vs.  SoC  data  estimated  by  Reynier  et  al.  as  a  close  approximation 
[  17  ].  The  ASC  vs.  SoC  and  ASa  vs.  SoC  sets  estimated  by  the  Yazami— 
Reynier  team  are  discrete  values  measured  at  large  intervals  of  SoC. 
We  fit  their  data  to  a  polynomial  approximation  (Fig.  6)  and 
interpolated  the  data  needed  in  the  required  intervals  of  SoC  for  the 
4.4-Ah  cell  discharging  at  4.4-A  rate. 

3.4.  Contributions  by  Rs,  Ra,  Rc.  ASC  and  ASa  to  heat  generation 

The  major  goal  of  this  work  has  been  to  demonstrate  an  approach 
to  measure  the  individual  contributions  by  the  five  modes,  Rs,  Ra,  Rc, 
ASC  and  ASa,  to  the  generation  of  thermal  energy  during  the 
discharge  of  a  Li-ion  cell.  Three  of  the  five  modes  (resistances  due  to 
the  electrolyte,  Rs,  anode,  Ra,  and  cathode,  Rc )  were  obtained  using 
impedance  technique;  see  Section  3.2  for  details.  Thermal  energy 
generations  due  to  the  Joule  heating  were  computed  as  i2Rs,  i2Ra  and 
i2Rc,  respectively  for  the  electrolyte,  anode  and  cathode;  i  represents 
the  dc  discharge  current,  which  was  fixed  at  4.4  A.  Thermal  energy 
generations  by  the  two  other  modes,  the  anode  and  cathode  en¬ 
tropies,  were  obtained  from  the  published  literature  data;  see  Sec¬ 
tion  3.3  for  details.  Computing  the  entropy-related  thermal  energy 
also  required  the  anode  and  cathode  temperatures  that  were 
measured  using  the  phase-shift  dependency  on  temperature;  see 
Section  3.1  for  details.  The  respective  heat  generations  were 
computed  as  (-TaASa)  and  (-TCASC). 

Fig.  7a  through  e  provides  a  summary  of  the  rate  of  heat 
generated  (in  units  of  Watt)  during  the  constant-current  discharge 
(4.4-A  constant  current)  by  each  of  the  five  modes.  Each  figure 
shows  six  different  graphs  corresponding  to  the  heat  data  for  one  of 
the  five  modes,  at  six  different  Tenv:  -10  °C,  0  °C,  10  °C,  20  °C,  30  °C 
and  40  °C. 
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The  total  rates  of  heat  generated  (in  units  of  Watt)  by  all  the 
five  modes  at  each  of  the  six  Tenv  are  shown  in  Fig.  7f.  For  each 
Tenv>  the  area  under  the  curve  provides  Qsum.  the  integral  of  the 
heat  generated  (in  units  of  Watt  hour)  by  all  the  five  modes  during 
the  entire  course  of  discharge  at  that  Tenv  In  other  words,  at  each 

Tenv,  Qsum  =  f  i2Rsdt  +  f  i2R3At  +  f  z2J?cdt+  [\-TM3)dt 
Jo  Jo  Jo  Jo 


+  /  (-TASc)dt;  where  t  is  the  discharge  time;  and  i  is  the  total 
Jo 

time  to  discharge  from  the  fully  charged  state  to  2.7  V  cut-off, 
during  constant-current  discharge.  Fig.  8  shows  the  Qsum  vs.  Tenv, 
as  well  as  the  integral  of  the  thermal  energy  due  to  each  of  the  five 
modes  of  heating. 

The  trend  in  Qsum  vs.  Tenv  implies  that  the  thermal  energy,  Q 
output  during  discharge  decreases  with  increase  in  the  environ¬ 
mental  temperature;  and  for  the  cell  model  tested  (Boston  Power 
Swing  4400),  the  Qat  20  °C  is  identical  to  the  heat  generated  by  the 
electrolyte  resistance,  Rs  at  20  °C.  At  lower  Tenv,  the  heat  generation 
is  higher,  while  the  least  amount  is  generated  at  higher  Tenv. 


Note  also  that  in  Fig.  7a-f,  the  time  to  discharge  to  cut-off  (2.7  V) 
is  shorter  at  lower  temperatures.  This  is  due  to  the  cell  reaching  the 
cut-off  voltage  sooner  at  lower  temperatures.  This  also  implies  that 
even  though  the  cut-off  voltage  is  fixed  at  2.7  V  during  discharge  at 
all  Tenv,  a  fraction  of  the  stored  chemical  energy  remains  un¬ 
discharged  in  the  cell,  and  the  magnitude  of  the  un-discharged 
energy  is  higher  at  lower  Tenv;  (see  also  Section  5). 

The  implications  of  the  observations  posted  in  Figs.  7a— f  and 
Fig.  8  are  elaborated  next. 

4.  Manifestations  of  Rs,  Ra,  f?c.  ASC  and  A Sa  on  the  thermal 
behavior  during  discharge 

It  is  obvious  from  the  data  in  Figs.  5  and  6  that  any  Tenv,  the  five 
internal  properties,  Rs,  Ra,  Rc,  A Sc  and  ASa,  vary  widely  over  SoC  and 
the  time-of-discharge.  Therefore,  during  discharge: 

•  The  rate  of  heat  generation  by  every  one  of  the  five  modes  will 
vary  with  time. 
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•  The  variations  do  not  follow  any  particular  behavior  or  law 
(linear,  exponential,  parabolic,  etc.). 

•  It  is  therefore  infeasible  to  generate  predictive  models  for  the 
thermal  behavior  of  a  Li-ion  cell  under  discharge. 

•  Continuous  monitoring  of  Rs,  Ra,  Rc.  Tc  and  Ta  during  the 
discharge  is  the  only  way  to  measure  the  heat  generated  by  each 
mode  and  each  component  (anode,  cathode  and  electrolyte)  of 
the  cell. 

Some  of  the  nuanced  aspects  of  the  non-linear  dependency  of 
the  five  modes  of  heat  generation  on  SoC  and  temperature  are 
described  next  with  reference  to  Joule  and  entropy  heating. 


Within  the  45%  and  90%  SoC  range,  Rs,  Ra  and  Rc  all  show  a 
minima  implying  that  the  three  resistances  will  generate  the 
least  amount  of  heat  within  this  range  of  the  SoC. 

Rs  dominates  over  Ra  and  Rc,  therefore  the  electrolyte  is 
responsible  for  most  of  the  resistive  heat  generation. 


4.2.  Entropy-assisted  heating  and  cooling 

Comparing  -TCASC  and  -TaASa  with  each  other,  (Fig.  7d  and  e): 

•  The  anode  always  cools  due  to  the  entropy  change  in  the  anode; 
such  a  behavior  of  the  anode  is  common  to  all  Li-ion  cells  that 
uses  carbon  anode. 

•  On  the  other  hand,  the  LMO-based  cathode  behavior  is  more 
complex,  whose  ASC  changes  in  the  positive  and  negative  di¬ 
rections  over  the  SoC  range  of  full  charge  to  full  discharge  [18], 
Thus,  ASC  tends  to  heat  the  cathode  in  certain  ranges  and  cool  it 


4.3.  Influence  of  foule  and  entropy  on  T„  and  Tc 

The  sum  ( i2Rc  +  (-TCASC))  influences,  at  least  partially,  the  Tc,  the 
temperature  of  the  cathode;  similarly,  (i2Ra  +  (-TaASa))  influences 
Ta.  However,  the  resistive  heating  by  the  electrolyte  has  stronger 
role  in  determining  the  temperature  of  the  anode  and  cathode  for 
most  part,  though  not  entirely.  The  dominance  of  Rs  over  Rc  and  Ra, 
and  the  proximity  of  the  electrolyte  to  the  anode  and  the  cathode, 
do  play  a  role  in  increasing  Ta  and  Tc  during  the  initial  stages  of 
discharge.  However,  the  entropy-based  thermal  properties,  -TCASC 
and  -TaASa  do  play  a  major  role  in  modulating  Ta  and  Tc.  Note  from 
Fig.  7,  the  heating  rate  due  to  i2Rs  is  almost  constant,  increasing  by 
less  than  10%  over  the  entire  course  of  discharge.  On  the  other  hand, 
the  rate  of  heating  due  to  all  the  other  modes  varies  more  widely, 
and  while  cooling  in  some  places.  The  modulations  in  Tc  can  be 
explained  by  the  fluctuations  in  (i2Rc  +  (-TCASC))  vs.  SoC;  and  those 
in  Ta  by  the  changes  in  (i2Ra  +  (-TaASa))  vs.  SoC. 

Fig.  9a— f  shows  the  sums  of  (i2Rc  +  (-TCASC))  and 
(i2Ra  +  (-TaASa))  through  the  course  of  discharge  at  all  the  Tenv;  also 
presented  in  each  figure  are  the  Joule  heating  by  the  electrolyte 
(i2Rs)  data  for  the  purpose  of  comparison.  Comparing  them  with  the 
Tc  data  in  Fig.  3  demonstrates  the  influence  of  (i2Rc  +  (-TCASC))  on 
Tc.  Almost  at  all  Tenv,  past  2000  s  of  the  C-rate  discharge  (SoC  <45%), 
Tc  begins  to  decrease  from  its  maximum  value.  Starting  at  the 
2000  s  mark,  (i2Rc  +  (-TCASC))  has  a  negative  slope  or  a  decreasing 
rate  of  heating,  followed  by  a  negative  sign  or  cooling  past  2500  s 
(SoC  <  30%).  By  comparison,  even  though  (i2Ra  +  (-TaASa))  turns 
negative  during  discharge,  due  its  relatively  smaller  amplitude,  it 
does  not  influence  Ta  until  renv  >  30  °C.  In  summary: 


4.1.  foule  Heating 

Comparing  Rs,  Ra  and  Rc  with  each  other,  the  electrolyte 
resistance: 

•  Rs  is  the  only  parameter  that  is  almost  invariant  with  Tenv, 
staying  within  the  range  of  0.02  Q— 0.024  O  across  all  SoC  and 
the  Tenv  range  (-10  °C  and  40  °C). 

•  Rc  varies  most,  from  0.0006  Q  (at  Tenv  =  40  °C)  to  0.05  fi  (at 
Tenv  =  -10  °C).  At  Tenv  <  0  °C,  the  cathode  resistance  Rc 
is  comparable  to  and  even  higher  than  the  electrolyte  resis¬ 
tance,  Rs. 

•  Ra  has  the  least  magnitude  of  the  three,  never  increasing  above 
0.006  Q. 


•  The  electrolyte  resistance  is  least  influenced  by  SoC  and  Tenv.  and 
it  dominates  the  resistive  heating  and  contributes  most  to  Tc  and 
Ta  over  much  of  the  discharge. 

•  The  entropy  change  in  the  cathode  influences  the  Tc  most  below 
45%  SoC. 

•  Below  45%  SoC,  the  entropy  change  in  cathode  is  negative 
causing  entropic  cooling;  however,  only  sufficiently  to  cause  a 
small  drop  in  the  Tc;  i2Rc  and  i2Rs  continues  to  dominate  over 
(-TCASC),  keeping  the  cathode  hot. 

5.  The  identity  of  stored  energy  at  all  Tem  of  discharge 

Three  sets  of  discharge  energy  data  are  needed  to  account  for 
the  manifestation  of  the  chemical  energy  stored  in  the  cell:  the 


Qx  (Wh) 


Fig.  7.  Heat  generated  during  discharge  by  the  individual  modes,  Rs,  Ra,  Rc,  A Sc  and  ASa  are  shown  in  (a)-(e);  and  the  sum  of  all  the  five  modes  in  (f).  Each  figure  has  six  graphs 
corresponding  to  measurements  made  at  six  different  Tenv. 


Fig.  8.  The  integral  of  the  heat  generated  due  to  the  individual  modes  ( J  i2Rsdt, 
J  i2Radt,  J  i2Rcdt,  J  (-TASa)dt  and  J  (-TASc)dt,  and  their  sum 
(Qsum  =  J*  i2Rsdt  +  i2Radt  +  J’  i2Rcdt  +  J\-T\Sa)dt  +  J\-ThSc)dt)  through 
the  course  of  4.4-A  rate  discharge,  from  full  charge  to  the  2.7  V  cut-off,  at  the  six 
different  Tenv.  The  dots  are  the  integrals  obtained  from  experimental  measurements, 
and  the  lines  are  polynomial  fits. 


electrical  energy,  thermal  energy,  and  the  energy  that  remained 
un-discharged  due  to  the  imposition  of  the  2.7  V  discharge  cut¬ 
off  voltage  during  the  discharge.  The  electrical  energy  Edi  is 
measured  in  terms  of  cell  voltage,  current  and  time  during 
discharge.  The  Qsum  data  in  Fig.  8  is  the  thermal  energy  at  each 
Tenv-  In  order  to  measure  the  un-discharged  energy  after  dis¬ 
charging  at  each  Tenv,  following  the  discharge,  we  allowed  the 
cell  to  equilibrate  with  the  room  temperature  (23.5  °C)  for  12  h, 
and  then  re-discharged  the  cell  at  a  slow  rate  of  1.1 -A  to  2.7  V 
cut-off  voltage.  The  electrical  energy  obtained  during  the  second 
discharge,  Ed2  was  assumed  to  be  the  entire  un-discharged  en¬ 
ergy  left  behind  during  the  first  discharge;  the  thermal  energy 
generated  during  the  second  discharge  was  assumed  negligible. 
The  results  are  shown  in  Fig.  10  that  include  the  electrical  energy 
(Edi )  and  the  thermal  energy,  QsUm  (obtained  during  the  first 
discharge  at  each  Tenv).  and  the  residual  electrical  energy,  Ed2 
(obtained  during  the  second  discharge  at  23.5  °C).  Also  shown  in 
the  figure  is  the  Edeliverable:  Edl  +  Qsum  +  Ed2.  Note  that  Edeliverable 
is  nearly  a  constant  at  14.4-Wh  through  all  the  Tenv.  which  sug¬ 
gests  at  14.4-Wh  is  the  maximum  amount  of  energy  that  the  cell 
would  ever  deliver. 
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Fig.  9.  Heat  generated  by  the  sum  of  resistive  heating  and  the  entropy  changes  in  the  anode  (i2Ra  +  (-raASa))  and  the  cathode  (i2Rc  +  (-TaASc))  during  discharge  at  different  Tenv  as 
indicated  in  (a)-(f).  Also  shown  are  the  heat  generated  by  the  resistive  heating  of  the  electrolyte  (i2 Rs)  at  those  Tem.  Note  that  the  heat  data  at  Tenv  =  30  °C  and  40  °C  are  similar,  and 
not  identical. 


6.  The  disconnect  between  Tsurf  and  the  dynamics  of  heat 
generation 

Most  Battery  Management  Systems  (BMS)  use  surface-mounted 
temperature  sensors  for  the  thermal  management  part.  The  data 
below  show  the  potential  disconnect  between  the  surface  tem¬ 
perature,  Tsurf  and  the  thermal  events  in  progress  inside  the  cell 
under  discharge. 

First,  the  surface  temperature  lags  behind  the  heat-generating 
events  inside  the  cell  by  tens  of  seconds.  Fig.  11  shows  the  Ta,  Tc 
and  Tsurf  for  the  4.4-Ah  cell  discharging  at  11 -A  rate  in  a  5  °C  Teriv. 
The  Ta  and  Tc  react  to  the  current  pulse  almost  instantaneously;  the 
15-s  lag  between  the  current  pulse  and  the  two  temperatures  is 
caused  by  the  data  collection  rates  for  Ta  and  Tc.  The  Tsurf  increased 
sluggishly,  not  catching  with  Ta  and  Tc  even  after  60  s.  Second, 
during  discharge,  there  is  virtually  little  correlation  between  Tsurf 
and  Tenv,  while  Tc  and  Ta  correlate  well  with  Tenv;  see  data  in  Fig.  12. 


The  data  in  the  two  figures  suggest  that  the  role  of  surface- 
mounted  temperature  sensors  might  need  reevaluation  to  ensure 
reliable  thermal  management  of  rechargeable  Li-ion  batteries  used 
in  a  wide  range  of  equipment  including  electric  vehicles,  ships, 
commercial  and  military  planes  and  satellites,  especially  when 
safety  is  utmost  essential.  In  other  words,  surface-mounted  tem¬ 
perature  sensors  should  be  replaced  by  internal  temperature  sen¬ 
sors  for  accuracy  and  reliability. 

7.  Conclusions 

The  paper  describes  a  newly  created  capability  to  estimate  in 
situ  the  heat  generation  from  the  perspectives  of  the  individual 
components  of  a  Li-ion  cell,  namely  the  anode,  cathode  and  elec¬ 
trolyte.  This  capability  is  quite  unlike  that  of  the  conventional  cal¬ 
orimeters  that  are  routinely  used  to  report  the  heat  generated  by 
the  entire  cell,  not  their  individual  components.  The  sensors  that 
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electrical  energy  only  at  high  Te„„.  At  lower  Tenv,  a  part  of  the  stored  energy  is  left  un 
discharged  in  the  cell.  At  all  Tenv,  the  thermal  energy  outputted  during  the  C-rat 
discharge  is  only  a  fraction  of  the  deliverable  energy. 


-10  0  10  20  30  40  50 

Tenv  (°C) 


Fig.  12.  The  graphs  demonstrate  the  absence  of  a  non-monotonic  association  between 
Tsurf  with  Tenv  that  may  pose  challenges  in  the  use  of  Tsurf  to  predict  the  cell's  internal 
temperatures  while  its  surface  is  actively  cooled.  By  contrast,  the  relationship  between 
Tc  and  Tenv,  and  Ta  and  Tenv  are  monotonic.  Moreover,  the  direct  measurements  of  Tc  and 
Ta  remove  the  need  to  predict  the  internal  temperatures.  The  TJ13*  represent  Tsurf,  Tc 
and  Ta  data  measured  at  an  arbitrary  time  during  the  4.4-A  discharge;  the  dots  are 
experimental  data  and  the  lines  are  polynomial  fits. 


have  empowered  this  new  ability  lend  themselves  for  everyday  use 
to  track  fire  safety  in  lithium  batteries,  and  they  are  quite  unlike  the 
conventional  calorimeters  that  are  relegated  to  laboratory  use  only. 

The  surface  temperature  (rsurf)  of  the  cell  lags  so  far  behind  its 
internal  temperatures  (Ta  and  Tc)  that  the  relevance  of  monitoring 
TSUrf  to  protect  a  cell  against  situations  like  thermal  runaway  is 
questionable. 

Our  work  confirms  that  the  rate  of  heat  generation,  Qby  a  Li-ion 
cell  under  discharge  varies  with  discharge  time,  t,  even  if  the 
discharge  current,  i  and  the  temperature  of  the  cell’s  environment, 
Tenv  are  constant.  In  other  words,  Q.  vs.  t  is  nonlinear,  even  if  i  vs.  t 
and  Tenv  vs.  t  are  steady.  The  nonlinearity  is  caused  by  nonlinear 
dependence  of  the  five  modes  responsible  for  the  heat  generation 
on  temperature  and  SoC.  The  five  modes  are  driven  by  five  different 
material  parameters  of  the  cell:  the  resistances  due  to  the  elec¬ 
trolyte,  anode  and  cathode,  and  the  entropy  changes  in  the  anode 
and  the  cathode.  Each  parameter  is  dependent,  to  varying  degrees, 


responses  by  Ta  and  Tc  are  much  faster.  Note  that  the  rate  of  data  collection  for  Ta  and  Tc 
was  in  15-s  interval;  Tsurf  was  collected  at  a  much  faster  rate. 


on  temperature  and  SoC.  Heat  generated  by  some  modes  dominates 
at  low  Tenv,  while  others,  at  high  Tenv  Most  importantly,  heat  gen¬ 
eration  by  each  mode  during  discharge,  behaves  exactly  as  pre¬ 
dicted  by  their  respective  property’s  dependence  on  Tenv  and  SoC. 
This  predictable  behavior  sets  up  an  excellent  opportunity  to 
monitor  heat  generation  during  discharge  and  safeguard  the  Li-ion 
cell  from  thermal  runaway. 

In  a  Li-ion  cell  under  normal  discharge,  Joule  heating  due  to  Rs, 
Rc  and  Ra  contribute  to  most  of  the  thermal  energy  released.  Of 
these  three,  Rs  generates  most  of  the  heat,  and  Ra  the  least;  and 
contribution  by  Rc  varies  with  Tenv:  at  low  Tenv,  the  heat  generated 
by  Rc  is  comparable  and  even  higher  than  that  by  Rs,  and  at  high 
Tenv,  it  is  small,  comparable  to  that  of  Ra,  and  less  than  3%  of  Rs. 

Since  variations  in  Rs,  Ra,  Rc,  A Sc  and  ASa  with  time  of  discharge 
are  nonlinear,  and  do  not  follow  any  particular  behavior  or  law  such 
as  linear,  exponential,  parabolic,  etc.,  continuous  monitoring  of  Rs, 
Ra,  Rc,  Tc  and  Ta  during  the  discharge  is  the  only  way  to  measure  the 
heat  generated  by  anode,  cathode  and  electrolyte.  Fast  tracking 
impedance  technique  is  ideally  suited  for  monitoring  and  safe 
guarding.  This  technique  measures  the  electrolyte,  anode  and 
cathode  resistances  once  every  2  min  [20],  and  anode  and  cathode 
temperatures  once  every  1  s  [13],  It  is  a  suitable  technique  for 
monitoring  heat  generation,  while  the  cell  is  under  charge  and 
discharge.  The  findings  described  in  this  work,  and  in 
Refs.  [10,20,21]  demonstrate  a  non-intrusive  online  technique  to 
safeguard  Li-ion  batteries  from  venting,  fire  and  thermal  runaway. 
This  novel  approach  of  tracking  the  five  modes  of  heat  separately 
from  each  other  opens  a  new  path  for  in  situ  study  of  thermal- 
runaway  in  Li-ion  cells. 
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Glossary 


C:  capacity  of  a  cell  in  units  of  ampere-hour  (Ah) 

C-rate:  amplitude  of  current  rated  at  the  full  capacity  of  the  cell 
Pdeiivembie-  maximum  deliverable  energy  (thermal  +  electrical)  from  a  fully  charge 
cell 

Ed2'-  un-discharged  energy  left  in  the  cell  when  the  cell  reaches  the  2.7  V  cut-off  limit 
f:  frequency  in  Hz 
i:  discharge  current 

Ra:  resistive  impedance  of  the  anode;  it  is  a  temperature-dependent  parameter 
Rc:  resistive  impedance  of  the  cathode;  it  is  dependent  on  temperature  and  SoC 
Rs:  resistive  impedance  of  the  electrolyte;  it  is  weakly  dependent  on  temperature 
and  SoC 

i2R„:  Joule  heating  in  anode  in  units  of  Watt  (W) 
i%:  Joule  heating  in  cathode  in  units  of  Watt  (W) 
i2Rs:  Joule  heating  in  electrolyte  in  units  of  Watt  (W) 

LMO:  lithium  manganese  oxide 

tj>:  phase  shift  between  the  applied  ac  current  and  resultant  ac  voltage 

AS:  change  in  entropy 

ASa:  change  in  entropy  of  the  anode 

/ISC:  change  in  entropy  of  the  cathode 

i:  total  time  to  discharge  from  the  fully  charged  state  to  2.7  V  cut-off,  during 
constant-current  discharge 

<2:  heat  generated  by  the  cell  during  charging  or  discharging 

Qas.’  heat  generated  by  change  in  entropy 

Qas„  ■'  heat  generated  by  change  in  anode  entropy 

Q^:  heat  generated  by  change  in  cathode  entropy 

Qsum.’  summed  value  of  the  integral  of  heat  generated  by  each  mode 
((Qsum  =  ^i2Rsdt  +  J*  i2Radt  +  J\2Rcdt  +  J\-TASf)dt  +  J\-ThSc)dt) ) 
SoC:  state  of  charge 
t:  time  of  discharge 
T:  temperature  in  K  or  °C 
Ta:  temperature  of  the  anode 
Tc:  temperature  of  the  cathode 

Tenv’  temperature  surrounding  the  outer  surface  of  the  cell,  typically  the  temperature 
of  the  environmental  chamber 

Tsurf-  temperature  at  the  outer  surface  of  the  cell  measured  by  a  K-type  thermocouple 


